Three-body forces: From cold atoms to nuclei 
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It is often assumed that few- and many-body systems can be accurately described by 
considering only pairwise two-body interactions of the constituents. We illustrate that 
three- and higher-body forces enter naturally in effective field theories and are especially 
prominent in strongly interacting quantum systems. We focus on three-body forces 
and discuss examples from atomic and nuclear physics. In particular, we highlight the 
importance and the challenges of three-nucleon forces for nuclear structure and reactions, 
including applications to astrophysics and fundamental symmetries. 
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I. INTRODUCTION 

In this colloquium, we discuss recent advances, chal- 
lenges, and perspectives of three-body forces in nuclear 
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physics and related areas. We start with a brief overview 
of the history of the subject. 

The simplest non-relativistic system in which three- 
body forces can appear is the three-body system. The 
study of the three-body problem has a long history in 
physics. The gravitational problem of the eart h-moon- 
sun s ystem was first considered by Newton ([Newton 
(|l687h ). t was a central topic in mathematical physics 
from the mid 1700's to the early 1900's. In gravity, only 
two-body interactions between point masses are present. 
However, three-body tidal forces arise if extended ob- 
jects, such as planets, are treated as point particles. 

Three-body forces also play an important role in quan- 
tum mechanics and the quantum many-body problem. 
If they are not already present at a fundamental level, 
three- and higher-body forces appear in effective theories 
or in practical calculations, where the degrees of freedom 
and the Hilbert space have to be restricted. Typically, 
there is a hierarchy of these forces and two-body forces 
provide the main contribution with three- and higher- 
body forces giving smaller and smaller corrections. 

A well known example of such a three-body force 
in nuc lear physics is the Fu j ita-M iyazawa three-nucleon 
force (jFuiita and Mivazawa Its main contribu- 

tion arises from the virtual excitation of a A(1232) res- 
onance in processes involving three nucleons interacting 
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FIG. 1 Three-nucleon force arising from virtual excitation of 
a A(1232) degree of freedom. Solid (dashed) lines indicate 
nucleons (pions). 



via pion exchanges illustrated in Fig.[T] In atomic, molec- 
ular and optical physics and quantum chemistry, non- 
additive forces arise if Born-Oppenheimer potentials are 
calculate d by integrating out the e lectronic degrees of 
freedom ([Kaplan and Novaro (1993)). In the theory of 
strong interactions, Quantum Chromodynamics (QCD), 
three-body forces arise already at a fundamental level 
as the three-gluon vertex induces an interaction between 
three quarks. 

Our discussion here is guided by effective field theory 
ideas. The concept of resolution plays a key role in this 
context. In a scattering experiment, a particle beam with 
de Broglie wavelength A can only probe structures at a 
scale R > A. Similarly, in a general process with typical 
momentum scale /x only physics at momenta p < fj, (or, 
equivalently, distances R > is resolved. Effective 

theories and the renormalization group provide a method 
to use this observation for quantitative calculations. The 
resolution scale in an effective theory is controlled by the 
momentum cutoff A. Physics at momentum scales larger 
than the cutoff is excluded from the effective theory and 
encoded in effective couplings, so-called low-energy con- 
stants. These constants and the relative size of two- and 
higher-body forces turn out to be resolution dependent. 
If one starts with two-body forces only at high resolution, 
many-body forces will appear naturally as the resolution 
scale is lowered. These induced many-body forces cap- 
ture the contributions of successive two-body interactions 
which are separated by a distance below the resolution 
scale. 

The natural size of many-body forces is determined 
by the underlying scales of the theory and the external 
momentum scale /x. Effective field theories provide a con- 
venient and systematic scheme to construct and estimate 
the size of many-body forces. It is one aim of this col- 
loquium to exemplify these features. For example, chi- 
ral effective field theories in the Weinberg scheme have a 
clear hierarchy of many-nucleon forces. The current state 
of the art is to include two- and three-nucleon forces. In 
the future, however, the inclusion of four-body forces may 



also be required to achieve the desired accuracy. 

The quantitative importance of three-body forces is 
well established in light nuclei. Therefore, they should 
contribute significantly in heavier nuclei as well. How- 
ever, their inclusion in many-body calculations is com- 
putationally challenging and has only become feasible 
in recent years. We discuss three-nucleon forces at dif- 
ferent resolution scales and show that their inclusion is 
mandatory for nuclear structure calculations. Whether 
this scheme breaks down for heavy nuclei beyond a cer- 
tain mass number is an open question, but at present 
there are no indications of such a breakdown. 

The outline of the colloquium is as follows: We start 
with a review of the theoretical framework for three-body 
forces including an illustration of their scheme depen- 
dence. In Section IIII1 we discuss the role of three-body 
forces in the universal regime of large scattering-length 
systems and give examples from nuclear and cold atom 
physics. This is followed by a discussion of three-nucleon 
forces in chiral effective field theory in Section IIVI The 
application of such forces to many-body systems and 
their relation to electroweak processes is presented in Sec- 
tions [V] and I VII Finally, we give an outlook and discuss 
future opportunities. 



II. THEORETICAL APPROACHES TO THREE-BODY 
FORCES: DEFINITIONS, DEPENDENCE ON SCHEME 
AND FRAMEWORK 

Effective field theory (EFT) provides a general ap- 
proach to understand the low-energy behavior of a phys- 
ical system. The underlying principle was concisely 
formulated in IWeinbergl ( 19791 ): The most general La- 
grangian consistent with all symmetries of the underly- 
ing interaction will generate the most general S-matrix 
consistent with these symmetries. If this idea is com- 
bined with a power counting scheme that specifies which 
terms are required at a desired accuracy, one obtains a 
predictive low-energy theory. The expansion is usually 
in powers of a low-momentum scale M\ OWl which can be 
the typical external momentum, over a high-momentum 
scale Mhigh- To illustrate this idea, consider a theory 
that is made of two particle species, a light and a heavy 



one with M\ ow <C Mm g h- We focus on soft processes 
in which the energies and momenta are of the order of 
the light particle mass (the so-called soft scale). Under 
these conditions, the short-distance physics related to the 
heavy-particle exchange cannot be resolved. However, 
it can be represented systematically by contact inter- 
actions between light particles. Consider heavy-particle 
exchange between the light ones at momentum transfer 
q 2 -C M^igh- The corresponding tree-level expression for 
the scattering amplitude is simply g 2 rj (M^ b — q 2 ) with 
g the heavy-light coupling constant. It can be expanded 



in powers of q 2 /M^ ih as: 
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This expansion can be represented in the EFT. At low 
momentum transfer q 2 , the effects of the pole from the 
heavy-particle exchange in Eq. ([T]) are captured by a 
series of local momentum-dependent interaction terms 
reproducing the expansion in Eq. ([T]) term by term. 
This idea is closely related to the multipole expansion in 
classical electr odynamics and the renormalization group 
dWilsonl (Il983h l. 

The interactions in EFTs are represented by operators 
Oi that are monomials in the quantum fields ip in the 
general interaction Lagrangian, 
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(2) 



These operators can contain an arbitrary number of 
derivatives and/or fields but must respect the symmetries 
of the underlying theory. The derivatives are converted 
to momenta and generate the momentum dependence ex- 
emplified in Eq. ([T]). The coupling constants gi can be 
ordered according to their importance at low energies 
from their scaling with q ~ Mi ow and M n i g h- Operators 
with a larger number of derivatives or fields are usually 
suppressed. This is the basis of the power counting of 
the EFT. 

An illustrative example is given by the Lagrangian 



N 



4nt = X>WW> 



(3) 



which contains momentum-independent two-, three-, ... 
up to TV-body contact interactions of a nonrelativis- 
tic field ip. In a natural theory without any fine 
tuning of parameters, the dimensionful coupling con- 
stants gi scale with powers of the high-momentum or 
breakdown scale M n i g h- Dimensional analysis requires 
that gi ~ (l/Mhigh) 3 * -5 , such that A-body interac- 
tions are suppressed by (Miow/Mhigh) 2 ^ -2 ^ compared 
to N — 1 successive two-bod y interactions (see, e.g., 
lHammer and Furnstahll f|2000h ^ . This type of scaling 
analysis is the basis of the suppression of many-body 
forces in the Weinberg scheme mentioned in the intro- 
duction. 

The values of the coupling constants gi are deter- 
mined completely by on-energy-shell information, up to 
a well-defined truncation error. The exact relation, how- 
ever, is not unique and depends on the renormalization 
scheme. In the construction of the most general La- 
grangian, many-body forces arise naturally. These many- 
body forces have to be determined from many-body data. 

A fundamental theorem of quantum field theory states 
that physical observables are independent of the choice 
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FIG. 2 Illustration of the field redefinition in Eq. ((SJ that 
trades an off-shell two-body interaction for a three-body in- 
teraction. 



of fiel ds in a Lagrangian ([Coleman et Ml (|l969h : lHaad 
(1958)). Consequently, they are invariant under redefini- 
tions of the fields in the effective Lagrangian. Off-shell 
amplitudes, however, change under field redefinitions and 
thus are not observable. In systems with more than two 
nucleons, one can trade off-shell, two-body interactions 
for many-body forces. 

In a quantum mechanics framework, unitary 
transformations provide an alternative formalism 
to field redefinitions and lead to the same result 
( Amghar and Desplanques ( 19951 ): IPolvzou and Glockld 
()199(i )). This explains how two-body interactions related 
by unitary transformatio ns can predict different bin ding 
energies for the triton ([Afnan and Serduke (<1973l) 1 if 
many-body forces are not consistently included. 

We use a simple EFT model, to illustrate how field 
redefinitions can be used to shift strength from off-shell 
two-body interactions to on-she ll three-body interactions 
([Hammer and Furnstahl (|200Cj )): 



(4) 

where T> = idt + V 2 /(2m) is the free Schrodinger op- 
erator. The model has a two-body contact interaction 
with coupling constant 172 and an off-shell two-body con- 
tact interaction with coupling 77 which we assume to be 
small. Now consider a field transformation 

i> — ► [1 + r){i>H)\ ip , ^ ^ [l + ^t/j)]^ . (5) 

Performing this transformation and keeping all terms of 
order 77 we obtain a new Lagrangian: 

£ = ^Vi> - g 2 {i>H) 2 - 4r7 52 (V> t V) 3 + 0(V 2 ) , (6) 

where the off-shell two-body interaction has been traded 
for a three-body interaction. This is illustrated in Fig. [5] 
Off-shell interactions always contribute together with 
many-body forces and only the sum of the two is mean- 
ingful. It is thus not possible to determine off-shell inter- 
actions from experiment. 

On the other hand, if new interactions are generated 
from a given two-body interaction by field redefinitions 
or unitary transformations and only the two-body part 
is retained, many-body observables will depend on the 
interaction ev en if they genera t e the same two-body ob- 



servables (see iFurnstahl et al\ (|2001l ) for an explicit il- 



lustration in the above model). For example, nuclear 
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FIG. 3 Illustration for the resolution dependence of two- and 
three-body interactions. 



matter binding curves will depend on the off-shell part 
of the two-b ody interaction gen erating so-called "Co- 



ester bands" (Coester et al. (1970)). In the 1970s, it was 



proposed that comparisons of nuclear matter calculations 
could help to determine the "correct" off-shell behavior 
of the two-body interaction. From a modern perspective, 
it is clear that this is not possible and the bands will dis- 
appear if the full transformed Hamilitonian is used. 

One might argue that it should be possible to find a 
suitable representation of the theory where three-body 
forces vanish or are very small. This strategy could be 
used to minimize the computational effort in many-body 
calculations. As discussed below, this is indeed possible 
for the universal EFT but only at leading order (LO). In 
theories with more complex operator structure and long- 
range interactions such as the chiral EFT, however, it 
is doomed to fail from the start. The various operators 
contribute differently to different observables and there 
is no optimal choice for removing the contributions of 
three-body forces at the same time for all observables. 

As illustrated above, the interaction strength may be 
shifted from two- to many-body forces. This can also be 
seen by changing the momentum cutoff A in the regu- 
lators used in explicit calculations. Once the couplings 
g s : of the effective Lagrangian, the low-energy constants 
(LECs), have been adjusted to selected data, predic- 
tions for other low-energy observables should be inde- 
pendent of the choice for A. Obviously, this adjustment 
depends on A implying that also the interaction strengths 
of two- and many-body interactions vary with the cutoff 
and are not unique. The idea is illustrated in Fig. [31 
where iterated two-body interactions at short-distance 
scales of ~ 1/A are not resolved. Note that in practice 
the 3NFs generated in this way cannot be disentangled 
from the 3NFs at an initial scale, which will also have 
short-ranged (and other) contributions. It is therefore 
model-dependent to distinguish such "generated" from 
"genuine" 3NFs as is done often in the literature and, 



therefore, we will not distinguish 3NF mechanisms in this 
colloquium. 

For two-body forces, this can be implemented by 
renormalization group (RG ) equations for the potential 



renormanzatio n group (KXj ) 
(iBoener et all (l2007bl [201ol) l 



dA 



V A (12) = F(V A ,A). 



(7) 



The RG equation describes the evolution with A of the 
matrix elements Va(12) of the potential in momentum 
space. The function F is defined such that the on-shell 
T-matrix is invariant under changes of the cutoff for mo- 
menta below A. This equation can be integrated from 
large A to lower cutoff scales. By construction, all two- 
body observables up to momenta of the order of the cut- 
off are invariant. Beyond this, e.g., for processes involv- 
ing external probes and more particles, observables will 
depend on the cutoff. Complete RG invariance is only 
achieved when many-body forces and many-body cur- 
rents are included. In principle, a similar RG equation 
for three- and higher-body interactions can be formu- 
lated. This has been realized i n pract ice with the similar- 
ity RG (SRG) (jBoener et all (|2007al) ) or by taking EFT 
three- nucleon forces (3NF s) as a general low-momentum 
basis ( Nogga et al. ( 2004h ). where the LECs are adjusted 
to few-nucleon data at the lower cutoffs. If the resolu- 
tion scale A is not too low, the contributions of many- 
body forces obtained in this way are of the size expected 
in EFT and small compared to the two-nucleon (NN) 
force contributions, but they are still quantitatively im- 
portant in state-of-the-art computations. The variation 
of the cutoff then enables one to estimate contributions 
of higher-body short-range interactions. This will be dis- 
cussed in more detail in Section IIVI 

While the RG evolution is already an interesting tool 
to estimate many-body contributions to specific observ- 
ables, it becomes even more valuable in many-body cal- 
culations (see Section [V} , where RG transformations to 
lower resolution lea d to greatly enhanced convergence 
(|Bogner et all (|201fj )). 



III. UNIVERSAL ASPECTS: FROM COLD ATOMS TO 
LOW-ENERGY REACTIONS AND HALO NUCLEI 

As discussed in the previous section, the short-distance 
properties of a physical system are not resolved in low- 
energy observables. If no massless particles are present, 
all interactions appear short-ranged at sufficiently low 
energy. It is then possible to formulate an EFT with 
contact interactions. 

Particularly interesting is the case of strong interac- 
tions characterized by a large scattering length a. Such 
systems are close to the unitary limit of infinite scattering 
length. It is obtained by taking the range of the interac- 
tion to zero while keeping a two-body bound state fixed 



at zero energyQ In this limit, the two-body scattering 
amplitude is scale invariant and saturates the unitarity 
bound. Formulated as a challenge to test many-body 
methods (Bcrtsch (1999)), this limit turned out to be rel- 
evant for a variety of systems. It is historically interest- 
ing to note that an approximation corresponding to the 
unitar y limit was already used in (jBeth and Uhlenbeck 
to calculate the second virial coefficient of a Fermi 
gas. Ultracold atomic gases can be tuned to the vicin- 
ity of the unitary limit using Feshbach resonances, while 
neutron matter is close to this limit through a fine tuning 
in nature. This gives rise to novel many-body phenom- 
en a, such as the BEC-BCS crossover in ultracold atoms 



(|Giorgini et all (120081 )) and the "perfect" liquid ob served 



in heavy- ion collisions (jSchafer and Teanev ( 200 

Here, we use the unitary limit as a starting point for 
an EFT expansion for strongly interacting quantum sys- 
tems with short-range interactions. This universal EFT 
is applicable to any system close to the unitary limit, 
i.e., any system with short-range interactions and large 
scattering lengths. Examples include halo states in nu- 
clear physics, ultracold atoms close to a Feshbach reso- 
nance, and hadronic molecules in particle physics. The 
breakdown scale Mhigh of this theory is set by the lowest 
energy degree-of-freedom not explicitly included in the 
theory. In nuclear and particle physics, this is typically 
given by one-pion exchange. In ultracold atoms, -Mhigh 
is determined by the van der Waals interaction, but the 
details depend on the system. The typical momentum 
scale of the theory is Mi ow ~ 1/a ~ k. For momenta 
k of the order of the breakdown scale Mhigh or above, 
the omitted short-range physics is resolved and has to be 
treated explicitly. 

The universal EFT exploits the appearance of a large 
scattering length, independent of the mechanism gen- 
erating it. Because the dependence of observables on 
the scattering length is explicit, it allows to unravel uni- 
versal phenomena driven by the large scatterin g length 
such as universal correlations of observables ( Phillips 



(|l968h . iTionl (Il975f )h the Efimov effect (lEfimovl Ij970'l 'l . 
and limit-c y cle p hysics ( Braaten and Hammerl (|2003 ). 



Mohr et oD ( 20061 )). For reviews of applications to the 



Braaten and Hammer 



physic s of ultracold atom s, see 
(120061) and IPlatterl (<2009h . The applic ations in nuclear 
and pa rticle physics were discussed in lEpelbaum et al. 



(|2009al ) and lHammer and Platterl (|2010t ) 

Three-body forces play an important role in the uni- 
versal EFT and we discuss their contribution in three- 
and higher-body systems in detail below. In the simplest 
case of spinless bosons, the leading-order Lagrangian can 



1 In real physical systems the strict unitary limit of zero-range in- 
teractions can, of course, not be reached. In low-energy observ- 
ables, however, the finite range R is not resolved and corrections 
are small (of the order of R/a or Rk). 




FIG. 4 The bubble diagrams with the contact interaction <?2 
contributing to the two-body scattering amplitude. 



be written as: 



(8) 



Extensions to more complicated systems are straight- 
foward. The terms proportional to gi and gs correspond 
to two- and three-body contact interactions. The dots 
represent higher-order terms suppressed by derivatives 
and/or more fields. 

The renormalized values of the coupling constants gi 
and c/3 are matched to observables in the two- and three- 
body system. In the two-body system, one typically takes 
the S-wave scattering length. The exact relation between 
the coupling gi and the scattering length depends on the 
rcnormalization scheme. Because of this matching pro- 
cedure, the EFT provides correlations between different 
observables based on the hierarchy of scales in the sys- 
tem. Given one set of observables, another set can be 
predicted to a certain accuracy. Depending on the ex- 
perimental situation, these correlations can be applied in 
different ways. 

Since the scattering length is large, a ~ l/M\ ow , the 
leading co ntact interact i on <?? has to be r esum med to 
all orders (<Kaplan et all (|l998l ): Ivan Kolckl (|l999h ). The 
two-body scattering amplitude is obtained by summing 
the bubble diagrams with the g2 interaction shown in 
Fig. |U This summation gives the exact solution of the 
Lippmann-Schwinger equation for the gi interaction and 
reproduces the leading term of the effective range expan- 
sion. Higher-order derivative interactions, which are not 
shown explicitly in Eq. (j5]), generate higher-order terms 
in the effective range expansion. Since these terms are 
set by Mhigh, their contribution at low energies is sup- 
pressed by powers of Mi ow /Mhigh and can be treated in 
perturbation theory. The first correction is given by the 
S-wave effective range, ro ~ 1/Afhigh- 

The S-wave scattering amplitude to next-to-leading or- 
der (NLO) then takes the form 



Ti{k) = 



1 



-1/a — ik 



1 - 



r Q k 2 /2 
-1/a — ik 



(9) 



where k is the relative momentum of the particles and 
the dots indicate corrections of order (Mi ow /A/hi g h) 2 for 
typical momenta k ~ Mi ow . If a is large and positive, Ti 
has a bound state pole at k = i/a. This corresponds to a 
two-body bound state (dimer) with binding energy Bi — 
l/(2/ia 2 ), where [i is the reduced mass of the particles. 
As a — > 00, this bound state approaches the two-body 
threshold. 
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FIG. 5 The integral equation for the boson-dimer scatter- 
ing amplitude. The single (double) line indicates the boson 
(dimer) propagator. 



The universal EFT shows its full strength in the two- 
body sector when external currents are considered. In 
contrast to other approaches, the coupling to currents 
is straightforward and current conservation is satisfied 
at each stage of the calculation. Gauge-invariant few- 
body contact terms are generated naturally by writing 
the most general effective Lagrangian. Applications to 
a variety of electroweak processes in the two-nucleq n 
sector have been carrie d out (see Beane et~al\ ( 200lh : 
iBedaaue and van Kolckl (l2002h for more details). Re- 
cently, these methods have also be en applied to neutron- 
rich s ystems and halo nuclei. In Hammer and Phillips! 
f)20Tlh . e.g., the electric properties of the one-neutron 
halo nucleus n Be were investigated. While this nu- 
cleus is nominally an 11-body system, the properties of 
its ground and first excited st ate can be described i n 
the framework of the h alo EFT ( Bedaque et all (2003a); 



Bertulani et al 



(|2002h ). This EFT exploits the small 
binding energy of these two states compared to the typ- 
ical energy scales of 10 Be (binding and excitation ener- 
gies). Thus, 11 Be can be treated as an effective two-body 
system of the 10 Be core and a neutron. A similar strategy 
was applied to calculate the ra d iative neutron capture 



on a 7 Li core ([Rupak and Higal ([20111 ); iFernando et al. 
4201 1JQ 



We now proceed to the three-body system where the 
term proportional to g 3 in Eq. (j8]) contributes. From 
naive dimensional analysis one would conclude that the 
53 term is of higher order (cf. discussion in Sec. [IT]). This 
is indeed the case for two-component fermions where the 
Pauli principle forbids three fermions to be close together 
in an S-wave. In general, however, naive dimensional 
analysis fails for large scattering length a. Again, we 
focus on the case of identical bosons which already con- 
tains the main features of the problem. The simplest 
three-body process to be considered is the scattering of 
a boson and a dimer. The integral equation for boson- 
dimer scattering is shown schematically in Fig. [5] For 



total orbital angular momentum L = 0, it takes the form: 

T 3 (k,p;E) = ^-M(k,p;E) + - [ dqq 2 T 3 (k, q; E) 

x M (q,P; E) 

-1/a + V 3 ? 2 /4 -mE-ie' 

where the inhomogeneous term reads 



M(k,p; E) 



+ kp + p 2 — mE 
2kp \k 2 — kp + p 2 — mE 



■In 



g(A) 
A 2 



(11) 

and a momentum cutoff A has been introduced to reg- 
ulate the integral equation. All other three-body ob- 
servables can be extracted from the amplitude T 3 taken 
in appropriate kinematics. In Eq. (p~0|) . H determines 
the strength of the three-body interaction 53(A) = 
— 4m (?2(A) 2 77 (A)/A 2 . The magnitude of the incom- 
ing (outgoing) relative momenta is k (p) and E — 
3k 2 /(4m) — l/(ma 2 ). The on-shell point corresponds to 
k = p and the scattering phase shift can be obtained via 
fccot<5 = l/T 3 (k,k; E)+ik. 

For H = and A -> 00, Eq. (JTO]) reduces to the STM 
equation first derived by Skorniak o y and Ter-Martirosian 
( Sko rniakov and T cr-Martirosia 3 |l957h which has no 
unique solution (Daniloy ( 196lh ) . The regularized equa- 
tion has a unique solution for any given (finite) value of 
the cutoff A but three-body observables show a strong 
dependence on the cutoff A. Cutoff independence of the 
amplitude is restored by an appropriat e "running" of 



H(A) w hich turns out to be a limit cycle (jBedaque et al. 
(|l999allbh ): 



H(A) 



cos[so ln(A/A») + arctanso] 
cos [so ln(A/A*) — arctanso] 



(12) 



where so ~ 1.00624 is a transcendental number and A* 
is a dimensionful three-body parameter generated by di- 
mensional transmutation. Adjusting A* to a single three- 
body observable allows to determine all other low-energy 
properties of the three-body system. Note that the choice 
of the three-body parameter A* is not unique and there 
are other definitions more dir ectly related to experiment 
(Bra aten and Hammer (Sxl)). 

The physics of this renormalization procedure is illus- 
trated in Fig.|6]where we show the unrenormalized three- 
body binding energies B3 in the case of positive scattering 
length as a function of the cutoff A (solid line). As the 
cutoff is increased, B 3 increases. At a certain cutoff (in- 
dicated by the dotted line), a new bound state appears 
at the boson-dimer threshold. This pattern repeats ev- 
ery time the cutoff increases by the discrete scaling factor 
exp(7r/so). Now assume that we adopt the renormaliza- 
tion condition that the shallowest state should have a 
constant energy given by the dashed line. At small val- 
ues of the cutoff, we need an attractive three-body force 
to increase the binding energy of the shallowest state as 
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FIG. 6 (Color online) Unrenormalized three-body energies 
Bz as a function of the momentum cutoff A (solid lines) . The 
dotted line indicates the cutoff where a new three-body state 
appears at the boson-dimer threshold (dash-dotted line). The 
dashed line shows a hypothetical renormalized energy. The in- 
set shows the running of the three-body force 33(A) ~ —H(A) 
with A. 



indicated by the arrow. As the cutoff is increased further, 
the required attractive contribution becomes smaller and 
around Aa = 1.1 a repulsive three-body force is required 
(downward arrow). Around Aa = 4.25, a new three- 
body state appears at threshold and we cannot satisfy 
the renormalization condition by keeping the first state 
at the required energy anymore. The number of bound 
states has changed and there is a new shallow state in 
the system. At this point the three-body force turns 
from repulsive to attractive to move the new state to the 
required energy. The corresponding running of the three- 
body force with the cutoff A is shown in the inset. After 
renormalization, the first state is still present as a deep 
state with large binding energy, but for threshold physics 
its presence can be ignored. This pattern goes on further 
and further as the cutoff is increased. 

The three-body force in Eq. (fl~2j) has exactly the right 
behavior to implement the strategy from the previous 
paragraph. Moreover, it breaks the scale invariance in 
the unitary limit, because the three-body parameter A* 
now provides a scale. However, due to the specific form 
of Eq. (fT2")l . a discrete scale invariance survives. Scaling 
transformations with the scaling factor Ao = exp(7r/so) 
leave H(A) and, consequently, three-body observables in- 
variant. This discrete scaling symm etry is the signature 
of an RG limit cycle ( Wilson ( 197ll )). In the three-body 
bound-state spectrum it becomes manifest through the 
Efimov effect: The appeara nce of a geome tric spectrum 
of three-body bound states (jEfimov 

The Efimov spectrum is illustrated in Fig. [Jj We show 
the energy variable K = sign(E)y/m\E\ as a function of 
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the inverse scattering length 1/a. The hashed areas indi- 
cate the three-atom (a < 0) and atom-dimer thresholds 



FIG. 7 (Color online) Illustration of the Efimov spectrum: 
The energy variable K = sgn(B) wm\E\ is shown as a func- 
tion of the inverse scattering length 1 /a. The solid lines indi- 
cate the Efimov states while the hashed areas give the scat- 
tering thresholds. The dashed vertical line indicates a system 
with fixed scattering length. 



(a > 0) where the Efimov states become unstable. The 
spectrum is invariant under the discrete scaling trans- 
formations K — > XoK and 1/a — > Xo/a. As a conse- 
quence, there is an accumulation of Efimov three-body 
states at the origin. The scaling symmetry relates Efi- 
mov states along any ray with fixed angle £ (cf. Fig. [7]). 
In general, these states correspond to different scattering 
lengths. A physical system with fixed scattering length 
is illustrated by the vertical dashed line. For fixed a, the 
discrete scaling symmetry is only manifest in the unitary 
limit 1/a = 0. 

The parameter A* can be used to set one of the three- 
body energies. All other states then follow from the dis- 
crete scaling symmetry. This explains why one param- 
eter is sufficient for renormalization of the whole spec- 
trum. The discrete scaling symmetry predicts infinitely- 
deep thr ee-body states. This is known as the Thomas 
collapse ( Thomas] (|l935l )V Physically relevant, however, 
are only states with energies \E\ <C M^ igh /m. All deeper 
states are ultraviolet artefacts of the effective theory and 
should be discarded. 

The discrete scale invariance also manifests itself in the 
log-periodic dependence of scattering observables on the 
scattering length. This scaling be havior has been con- 
firmed in cold atom experiments ([Ferlaino and Grimm 
(|20inh ). ja such experiments, the scattering length can 
be varied using Feshbach resonances. The scattering- 
length dependence of three-body recombination rates 
provides indirect information on the Efimov spectrum. 
For negative scattering length, the Efimov states hit 
the three-atom threshold, E = 0, for certain values 
of a (cf. Fig. [7]) and lead to enhanced recombination 
rates. For positive scattering length, the Efimov states 
become unstable already at the atom-dimer threshold, 
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FIG. 8 (Color online) Three-body loss coefficient Kz in a gas of ultracold 7 Li atoms as a function of scattering length (in units 
of the Bohr radius ao) for the \uif = 1) state (red solid circles) and the \m,F = 0) state (blue open diamonds). The solid lines 
represent fits to th e universal EFT pr ediction. The dashed lines represent the K3 ~ a 4 upper (lower) limit for a > (a < 0). 
(Figure taken from lGross et al\ ||201Ct ).1 



E = — I /(ma 2 ), but interference effects lead to minima 
and maxima in the rate at E = 0. Ideally one would like 
to see multiple recombination features on each side of the 
Feshbach resonance. For equal mass particles, this is not 
a simple task because of the large scaling factor. When 
effective-range effects are included perturbatively as in 
Eq. ([9]) , the discrete scale invariance is softly broken, but 
the effects of t he break i ng on the recombination rate can 
be calculated (jji et all l|2012h ). 

As an example, we show in Fig. [5] the three-body loss 
coefficient K3 in a gas of ultracold 7 L i atom s measured 
by the Khaykovich group (|Gross et all (|2010h ) as a func- 
tion of scattering length (in units of Bohr radius ao) for 
the \mp = 1) state (red solid circles) and the |m_F = 0) 
state (blue open diamonds) . The data show that the po- 
sitions and widths of recombination minima and Efimov 
resonances are identical for both states, which indicates 
that the short-range physics is nuclear-spin independent. 
The solid lines giv e fits to the analytical expre ssions of 
the universal EFT ( Braaten and Hammer ( 20061 )) and re- 
produce the data very well. 

A more direct way to observe Efimov states is to 
populate these states directly through radio frequency 
transitions. This is difficult because of their short 
lifetime and has on l y rec e ntly been achie v ed fo r 6 Li 
atoms ( Lompe et al. ( 2010l ): Nakaiima et al. ( 201ll )). 

The integral equations for the three-nucleon problem 
are a generalization of Eg. (fTTJl). (Fo r their explicit 
form and derivation, see iBedaque et all (|2003bJ l.) The 
leading-order three-body force is required in all chan- 
nels where short distances are not shielded by the an- 
gular momentum barrier and/or the Pauli principle. For 
S-wave nucleon-deuteron scattering in the spin-quartet 
channel the three-body force is of higher order, and the 



spin-quartet scattering phases can therefo re be predicted 



to hig h precision from two-body data (jBedaque et al. 
(|l998l ): IBedaque and van Kolckl (|l998r i). In the spin- 



doublet channel there are two coupled channels but the 
renormalization is similar to the three identical-boson 
case. Thus, one needs a new parameter which is not 
determined in the NN system in order to determine the 
(leading) low-energy behavior of the three-nucleon sys- 
tem in this channel. A comprehensive discussion of three- 
bo dy_Jo£ce_effects in the three-nucleon system was given 
bv lGriefihammerl (|2005| ). 

The three-body parameter gives a natural explanation 
of universal correlations between different three-body ob- 
servables such as the Phillips line: a correlation between 
the triton binding energy a nd the sp i n-dou blet neutron- 
deuteron scattering length (j Phillip 1 dl968h ). These ob- 
servables are calculated for different two-body potentials 
that reproduce the NN scattering phase shifts but the 
three-body parameter is not constrained by the data. 
This generates a one-parameter correlation between dif- 
ferent three-body observables. These correlations are 
driven by the large scattering length and are indepen- 
dent of the mechanism responsible for it. As a con- 
sequence, they occur in atomic syst e ms su ch as 4 He 
atoms as well (jBraaten and Hammer (2006) ). Fo r an 
overview of this topi c see Epelbaum et all (|2009al ) and 



Hammer and Platterl (2010) 



The universal EFT has also been applied in the four- 
body sector. A study of the cutoff dependence of 
the four-body binding energies revealed that no four- 
body fo£ceJs_j^guired_for renor malization at leading or- 
der (jPlatter et all (|2004 l20Q5j) ). Thus, the four-body 
force is a higher-order effect. As a consequence, there 
are also universal correlations in the four-body sec- 
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four-body parameter for shallow four-body states. Their 
prediction of the resonance positions lead to the experi- 
mental observ ation of univer s al tet ramer states in ultra- 
cold caesium ( Ferlaino et al. ( 20091 )). This, in turn, has 



B( H) [MeV] 



FIG. 9 (Color online) The Tjon line correlation between 
B( 3 H) and _B( 4 He). The experimental value is shown by 
the cross. Th e shad ed band gives the LO universal result 
(jPlatter et al\ (|2005l )). while the outer solid lines include 
NLO correction s obta ined using the resonating group method 
(jKirscher et all (|2010l )). The pluses and diamonds show cal- 
culations using p henomenological N N and NN+3N poten- 
tials, respectively ijNogga et all (120001^ 'I : the circle gives a chi- 
ral EFT result at N 2 LO (|Epelbaum etaH fc00$ )); and the 
triangles are based o n an SRG-evo l ved N 3 LO NN potential 
(jNogga et all l|2004 ); lHebeler et all <|201ll )1. 



tor driven by the large sca ttering length . The prime 
example is the Tjon line (|Tionl |l975)): a correla- 
tion between the triton and alpha-particle binding en- 
ergies, B( 3 H) and f?( 4 He). Higher-order range correc- 
tions break the correlation and generate a band. In 
Fig.[9l we show this band together with calculations usin g 
phenomenological NN potentials (jNogga et al. 
a chiral NN potential at next-to -next-to- leading order 
(N 2 LO) (lEpelbaum et all (|2002h . SRG-evolved next- 
to-next-to-next-to-leading order (N 3 LO) NN potentials 
(|Nogga et al\ (|2004l k iHebeler et all (|201ll )). and the ex- 
perimental value. All calculations with interactions that 
give a large scattering length must lie within the band. 
Different short-distance physics and/or cutoff depen- 
dence should only move the results along the band. This 
can be observed explicitly in the results for the SRG- 
evolved N 3 LO NN potential indicated by the triangles 
in Fig. as well as in few-bo dy calculations with low- 
momentum interactions Mowfe ( Nogga et all ( 2004 )). 

The absence of a four-body force at leading or- 
der also implies a univers al four-body spectrum. In 
lHammer and Platter! ((2007h the dependence of the four- 
body bound-state spectrum on the two-body scattering 
length was investigated in detail and summarized in a 
generalized Efimov plot for the four-body spectrum. In 
particular, it was found that there are two four-body 
states tied to every Efimo v trimer. In a subsequent study, 
von Stecher et al. (|2009h extended these calculations to 



led to increased theoretical activity in this area. The sen- 
sitivity of tetramer energ ies to a four-body scale was, for 
example, investigated bv lHadizadeh et all (|201l! ). Four- 
body recombination and o ther scattering processes were 
calculated by Deltuva (see lDeltuval <|2012l ) and references 
therein). 

The bound-state properties of larger systems of bosons 
interacti ng through short-range in teractions were consid- 
ered bv lHanna and Blumd ( 2006! ). Using Monte Carlo 
methods they showed that universal correlations between 
binding energies can also be obtained. Calculations for 
larger number of particles using a model that incorpo- 
rates the universal behavior of t he th ree-body system 
were carried out by von Stecher ( 2010l ). These findings 
indicate that there is at least one iV-body state tied 
to each Efimov trimer and numerical evidence was also 
found for a second excited 5-body state. In a subsequent 
study (jvon Steelier! (|201l|)), the energies and structural 
properties of bosonic cluster states up to N — 6 were 
calculated for various two-body potentials. Besides the 
lowest cluster states, which behave as bosonic droplets, 
cluster states bound weakly to one or two atoms forming 
effective cluster-atom "dimers" and cluster-atom-atom 
"trimers" were identified. Thus the prospects for ob- 
serving universal physics in larger few-body systems are 
excellent. 

Recently, a geometric spectrum of universal three-body 
states has also been predicted for atoms with dipolar in- 



teractions ( Wang et al 



the four-particle threshold and confirmed the absence of a 



In this case, the struc- 
ture of the interaction is very similar to the nuclear tensor 
force generated by one-pion exchange. If the dipole mo- 
ments of the atoms are aligned, the interaction is attrac- 
tive in a head-to-tail configuration of the atoms and re- 
pulsive side-by-side, like for dipole magnets. If the dipole 
moments are anti-aligned, the interaction is opposite, re- 
pulsive and attractive, respectively. This might open the 
possibility to simulate the nuclear tensor force in exper- 
iments with ultracold atoms. 



IV. THREE-BODY FORCES IN FEW-NUCLEON 
SYSTEMS 

Three-body forces are especially important in nuclear 
physics. Phenomenological studies indicate, e.g., that 
the contribution of 3NFs to binding energies of light nu- 
clei is quantitatively signi ficant, of the order 20% (see 
Pieper and Wiringal (|200lh ). 

In nuclear systems, the long-range parts of nuclear 
forces are mediated by pion exchanges in addition to 
short-range contact interactions discussed in the preced- 
ing section. The pion mass is comparable to the momenta 
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FIG. 10 Topology of the two-pion-exchange 3NF. Solid 
(dashed) lines indicate nucleons (pions). 




FIG. 11 Topology of the leading mid-range (left) and short- 
range (right) 3NFs. Solid (dashed) lines indicate nucleons 
(pions) . 



in typical nuclei. Therefore, it cannot be expected that 
pionless EFT is applicable, unless one considers very spe- 
cific observables, e.g., for halo nuclei where the nucleon 
separation energy presents a low scale. 

An important signature of 3NFs is the model depen- 
dence of few-nucleon predictions when only NN interac- 
tions are employed in the calculation. This is demon- 
strated by the 3 H and 4 He binding energies in Fig. |9] 
The model dependence indicates that the missing three- 
nucleon and higher-body int eractions are different for 
each NN potential employed (lAmghar and Desplanquesl 
(|l995h : IPolvzou and Glockld jl990j)). Fortunately, 4N 
and higher-body forces are expected to be further sup- 
pressed. We will come back to this issue, but assume for 
the moment that such contributions are negligible. 

It is therefore required to formulate both two- and 
three-body forces within one systematic scheme. Histori- 
cally, this has not been the case. In most models, the 
main contribution is related to t he two-pion-exchange 
contri bution depicted in Fig. [TUJ iFuiita and Mivazawal 
(|l957l) realized that this model can be constrained using 
pion-nucleon scattering data and found that the interac- 
tion is dominated by P-wave pion-nucleon (7rN) interac- 
tions. This was the birth of modern 3NFs which were 
most ly developed indep endently of NN interactions (see, 
e.g., ICoon et all dl979h l. Two-pion- e xchan ge 3NFs are 
generally of the form (see lFriar et al. ( 19991 )): 
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where k label particles and a,/3,7 isospin, is the 
pion mass, f v — 92.4 MeV the pion decay constant, and 
gA the axial pion-nucleon coupling. The constants Cj 
are different for all models. In these equations, we have 
neglected cutoff functions that are required to regularize 
3NFs at short distances. 



In these models, the 3NFs are unrelated to the NN in- 
teraction, which shows up in a strong model dependence 
of predictions based on combining such 3NFs and dif- 
ferent NN interactions. Although often parts of the pa- 
rameters are adjusted using 7rN scattering data, one still 
needs to adjust a parameter of the 3NF, e.g., a cutoff 
parameter, such that the prediction for the 3 H binding 
energy agrees with experiment. Such combinations are 
not based on a consistent framework. They do not de- 
scribe all available 3N scattering data, but they improve 
the de scription of many low-energ y few- nucleon observ- 



ables ( Kalantar-Navestanaki et al. ( 2012h ). In part, this 



is related to the universal correlations of observables dis- 
cussed in the last section. Such observables are thus not 
useful to pin down the spin-isospin structure of 3NFs. 

Therefore, investigations have concentrated on 
intermediate-energy nucleon-deuteron scattering. Us- 
ing phenomenological forces, it can be shown that, 
for this energy range, observables exi st that are 
sensit ive to the structure of 3NFs (jWitaia et al. 
(|200ll) ). Due to a series of sophisticated nucleon- 
deuteron scattering experim e nts, d ata is now available 



( Kalantar-Navestanaki et al. (|2012l )). Unfortunately, 



this data has not been analyzed yet in a framework that 
provides consistent NN and 3N interactions. Comparison 
of the data to the predictions based on phenomenological 
forces show that the current models do not describe the 
intermediate-energy data very well. We note that im- 
provements of the models ha v e been suggested (see, e.g., 



Kalantar-Navestanaki et al. ( 2012t ) for more details) 



Here, we will focus on selected low-energy observables 
within a systematic approach to nuclear forces. 

Such a systematic approach has bee n developed base d 
on chiral EFT (for recen t reviews, see lEpelbauml (120061) : 
Epelbaum eTal\ (|2009al ): iMachleidt and Enteml (|201ll )). 
Based on the symmetries of the QCD Lagrangian, most 
importantly its spontaneously broken chiral symmetry, 
it is possible to formulate an EFT in terms of nucleon 
degrees of freedom and the nearly massless Goldstone 
bosons of QCD: the pions. Symmetry considerations suf- 
ficiently constrain the interactions of pions with them- 
selves and with nucleons to develop a systematic power 
counting scheme: chiral perturbation theory (ChPT). 
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The expansion parameter is Q/A x , where Q ~ is 
a typical momentum or the pion mass and A x w 1 GeV 
is the chiral-symmetry breaking scale. The terms of the 
Lagrangian that are relevant to our discussion here read 



sion, two other topologies (see Fig. [TT]) contribute 
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(15) 



where ip and 7? are nucleon and pion fields, respectively. 
The Lagrangian includes as 7T7tNN vertices the same 
LECs ct of the two-pion-exchange 3NFs. In addition, 
these LECs also contribute to the subleading two-pion- 
exchange NN interaction, which shows the strong con- 
nection of NN and 3N forces. 

The challenge for nuclear forces is that, because of 
bound state, parts of the interactions are nonpertur- 
bativc in contrast to the interactions in pionic or 7rN 
systems. This issue has been tackled assuming that 
the power counting can be applied to all non-reducible 
diagrams without purely nucleonic intermediate states. 
These diagrams form a potential that is the n summed 
to all orders solving the Schrodinger equation ( Weinberg] 
l|l99Ch ). The approach requires numerical regulariza- 
tion of the potential introducing a regulator dependence. 
In state-of-the-art applications, cutoffs are presently re- 
stricted to values A < 450 - 600 MeV. There is a 
lively discussion whether this constraint is an artifact 
and can be overcome by an improved power counting, or 
whether it is an inherent requirement of a nonperturba - 
tive extension o f ChP T ( Epelbaum and Gegelia ([2009); 
Hammer et all (|2007t) ; iNogga et all (|2005l) ; IValderrama 
(|201ll) V For applications to many-nucleon systems, lower 
cutoffs are advantageous and in some cases, computations 
are only feasible with lower cutoffs. For the estimate of 
missing higher-order contributions, we will assume a high 
scale A x « A in this section. 

In chiral EFT, NN, 3N and higher-body forces can be 
derived consistently. The general result is that higher- 
body forces are suppressed compared to lower-body ones. 
This justifies our assumption that 4N and higher-body 
forces are further suppressed. For NN interactions, one 
finds that the longest-range part is one-pion exchange, 
which is also the basis of state-of-the-art NN models. 
The first 3NF co ntribution is suppressed by (Q/A x ) 3 



([van Kolckl (|1994l )) and contains the two-pion-exchange 
part given by Eq. (ITS")) . At the same order of the expan- 
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Usually, these parts are called the D- and i?-term. The 
D-term is of mid range (one-pion-exchange-short-range) 
and the E-term is of short range. This implies that the 
.E-term coupling can only be obtained from few-nucleon 
observables, whereas the D-term strength is also related 
to weak or pionic processes involving two nucleons (see 
Section |VI|) . Following standard conventions, we intro- 
duce two dimensionless couplings c D = D f%A x and 
ce = E f*A x . As noted above, due to the a vertices 
of the Lagrangian, ChPT provides relations between the 
strength of the two-pion-exchange NN interaction and 
ySNF _ Thjg i eve i f consistency can only be implemented 
in the framework of ChPT. For the results given here, 
the A is not treated as an explicit degree of freedom. 
Since the mass difference of the nucleon and the A is only 
~ 2m w , an explicit inclusion is exp ected to improve th e 
convergence of the c hiral expansion (IKaiser et al. l (ll998T ): 



Krebs et all (|2007t ): lOrdonez et all (| 19941 ) ). ~Fbr 3NFs, 



3NF 



the leading A contribution is entirely included in V 2 
and shows up in larger strength constants C3 and C4 en- 
hancing the two-pion-excha nge contribution s comp ared 
to the other two topologies (jEpelbaum et all (|2008l )). 

In nuclear systems the separation of the high and the 
low scales (given by A and the pion mass or a typi- 
cal momentum of the system, respectively) is not ex- 
ceedingly large, which implies a slowly converging chiral 
expansion. Especially for intermediate-energy nucleon- 
deuteron scattering, the expansion parameter is esti- 
mated to be ~ 1/2 or larger. Therefore, calculations up 
to order Q 3 (including the leading 3NFs) are useful only 
up to nucleon laboratory energies of ~ 100 MeV. Fortu- 
nately, t he Q 4 3NF con t ributions ha ve been completed 



201l[ )L and applications 



recently ([Bernard et all ([2008, 
are under way. 

Before calculations based on chiral 3NFs can be per- 
formed, one needs to determine the LECs c,, cd and ce- 
The Ci constants have been determined from NN data as 
well as 7rN data. The results are summarized in Table UJ 
For simplicity, we have omitted the theoretical uncer- 
tainties and only give the central values. Most determi- 
nations are in agreement within the uncertainties, but 
deviations of the different determinations can be sizable, 
of the order of 30%. For our purpose here, this accuracy 
is sufficient and comparable to higher-order contributions 
that we do not take into account. This problem will be- 
come more relevant, when the subleading parts of the 
3NF will allow us to increase the accuracy of our predic- 
tions. In principle, these constants can be obtained in- 
dependently of the NN interaction. So their size should 
not depend on the regulator chosen or on the specific 
realization of chiral NN potentials. 
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TABLE I Comparison of different d determinations. The a's 
are given in GeV" 1 . At present, the determinations using NN 
observables require further constraints from 7rN observables to 
be conclusive. The last column indicates whether the c, values 
are mostly based on NN or 7rN data. We als o show the re- 
sults b ased on resonance saturation (res) from lBernard et al\ 
(|1997| ) (note that we omitted the 7rN fit from that paper). 



TABLE II Power counting predictions and explicit results for 
the binding energy B and the expectation values of NN and 
3N forces for 4 He. Cutoffs and energies are given in MeV. 





Cl 


C3 


C4 




Fettes et al. (1998) (Fit 1) 


-1.2 


-5.9 


3.5 


7rN 


Biittiker and Meifiner (2000) 


-0.8 


-4.7 


3.4 


7rN 


Meifiner (2007) 


-0.9 


-4.7 


3.5 


7rN 


Rentmeester et al. (2003) 


-0.8 


-4.8 


4.0 


NN 


Entem and Machleidt (2002) 


-0.8 


-3.4 


3.4 


NN 


Entem and Machleidt (2003) 


-0.8 


-3.2 


5.4 


NN 


EDelbaum et al. (2005) 


-0.8 


-3.4 


3.4 


NN 


Bernard et al. (1997) 


-0.9 


-5.3 


3.7 


res 



Based on naturalness arguments, one would expect 
that the q's are of the order of A^ 1 ~ 1 Gev -1 . It sticks 
out that C3 and C4 are larger than this estimate. This can 
be understood based on resonance saturation, where the 
large Cj's are related to the small A to nucleon mass d if- 
ference ~ 1/(toa — m) (see, e.g., Bernard et al. ( 19971 )). 
Taking A's explicitly into account reduces the magnitude 
of the Cj considerably so that an improved convergence 
of the chiral expansion can then be expected (see Sec- 
tion EED and iKrebs et all (|2007l ) in the context of NN 
interactions). 

Finally, we need to determine the constants cd and 
ce- Usually, combinations of ce and ce are found that 
make sure that the 3 H binding energy is described cor- 
rectly. Then different strategies have been used to con- 
strain ce from few-nucleon data, e.g., b y fitting the dou- 
blet n eutron-deuteron scattering le ngth ( Epelbaum et all 
(|2002f) ). the binding energy of 4 He (iNogga et all (I2006T)). 
or the radius of 4 He ( Navratil et aLl i 20071 )). In addition, 
the 3 H beta decay half-life can be used to constrain cd 
(see Section |VT|) . In particular, the fit of ce to the 3 H 
beta decay half-life or to the radius of 4 He have been 
shown to lead to a good overall description of light and 
p-shell nuclei. It is important to note that many low- 
energy observables are already well described once ce 
and ce combinations have been chosen that describe the 
3 H binding energy correctly. Therefore, the sensitivity of 
these observables on ce is low and a considerable uncer- 
tainty remains. Possibly, for higher-order calculations, 
other strategies need to be devised to obtain more accu- 
rate determinations of ce and eg. 

Since the separation of scales is not very large and since 
there are ongoing discussions on the size of the high scale 
for nuclei, it is instructive to calculate the contributions 
of NN and 3N forces to the binding energy of light nu- 
clei. These contributions are not observables, neverthe- 



A 


B 


(Vnn) 


(V3NF} 


1 <V3NF> 1 IXKl 

| (V NN > | L /0 1 


450 
600 


27.65 
28.57 


-84.56 
-93.73 


-1.11 

-6.83 


1.3 
7.2 



less their relative size can be estimated and compared 
to the power counting estimate. For this estimate, we 
use t he realization of c hiral EFT interactions at order 
Q 3 of [Epelbaum et all (|2005f ) (with A = 700 MeV). In 
this work, the NN potential has been fitted for different 
cutoffs, which can be used to investigate the scale depen- 
dence of chiral 3NFs. For the 3NF, we use the same Ci val- 
ues as for the NN part. The ce and ce values have been 
determined by a fit to the 3 H binding energy and the dou- 
blet neutron-deuteron scattering length. The chiral 3NFs 
have been regularized using a cutoff function depending 
on the relative momenta in the in- an d outgoing state 



where the cutoff is identical to the A of lEpelbaum et al 
(2005). The results are given in Table [TT1 As one can 



see, the binding energy of 4 He is close to the experimen- 
tal value of 28.30 MeV. The remaining deviation from 
experiment is comparable to the cutoff dependence and 
indicates the contribution that can be expected from or- 
der Q 4 . The leading 3NF is a Q 3 contribution. Assum- 
ing a typical momentum ~ m OT and A x = 500 MeV, we 
expect a contribution of approximately 2% to the po- 
tential energy. It is apparent that the contribution of 
3NFs strongly depends on the cutoff. For the first case 
in Table [Til the size is smaller than expected, which is 
no contradiction to the power counting. For the sec- 
ond case, the 3NF contribution is somewhat larger than 
naively expected. The estimate is still within a factor of 
3 — 4 correct (a natural-sized number), but it shows the 
enhancement of the 3NF due to the A resonance. 

In summary, the overall size of 3NF contributions is 
as expected from the power counting once the contri- 
bution of the A resonance has been taken into account. 
The deviation of the binding energy for 4 He can also be 
expected from a higher-order contribution. On a quan- 
titative level, this deviation indicates that high precision 
can only be expected for a Q 4 calculation. 

At order Q 3 , there are also nucleon-deuteron scatter- 
ing calculations available. At intermediate energies, the 
results are strongly dependent on the cutoff. For low en- 
ergy, however, many observables can be accurately pre- 
dicted. In the left panel of Fig. [TU we show as an ex- 
ample the elastic nucleon-deuteron cross section. For the 
elastic cross section, data and prediction are in excellent 
agreement, and the order Q 2 and Q 3 results are similar 
indicating that the calculation is converged with respect 
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FIG. 12 (Color online) Left panel: Elastic nucleon-deuteron cross section at 10 MeV. The almost indistinguis h able bands cor- 
respond to chiral Q 2 ( red, d a rk grey) and Q 3 (cyan , light grey) calculations. Data are from lHowell et all (|l987l ); lRauprich et al\ 
(|1988| ); TSagara et al\ (|1994T ); fSperisen et al\ (I1984T) . Righ t panel: Nucleon-deuteron breakup cross section at 19 MeV for the 
space-star configuration at a = 56° (see Lev et all (|2006T l for the definition of the kinematics). The bands are the same as in 
the left panel. The data is from lLev et all ( |2006T ). 
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FIG. 13 (Color online) The nucleon vector analyzing power 
A y for elastic nucleon-deuteron scattering at 10 MeV. The 
cyan band covers predictions based on chiral NN and 3N forces 
at order Q 3 for various cutoffs. The solid, purple line is the 
result of the CD-Bonn NN combined with the TM99 3N po- 
tential. Data as in the left panel of Fig. 1121 Repri nted with 
permission from iKalantar-Navestanaki et al\ (|2012T ) . 



to the chiral expansion. Whereas the bulk of the observ- 
ables at low energies are nicely reproduced, there are a 
some exceptions. One of them is a specific breakup con- 
figuration shown in the right panel of Fig. 1121 Again, the 
Q 2 and Q 3 results nicely agree indicating convergence of 
the chiral expansion and, therefore, small contributions 
of 3NFs. Unfortunately, there is a large discrepancy to 
the data. This is still an unresolved problem. 

A different example is the analyzing power at low en- 
ergy shown in Fig. [T3l The solid line represents the 
result based on high-precision phenomenological forces, 
which clearly disagrees with the data. This is a com- 
mon feature of all available calculations based on phe- 



nomenological forces a nd has been discussed vivi dly in 
the literature (see, e.g., Miller and Schwenk ( 2007|)'). A t 



(2005). 



order Q 3 , for the realization of lEpelbaum et al. 
the result seems to agree with the data. However, the 
cutoff dependence is unusually large for these small en- 
ergies and a detailed analysis reveals that this agreement 
in the three-body sector can be traced back to deficien- 
cies in the description of NN data. Therefore, at order 
Q 3 , the analyzing power cannot be properly predicted 
but is merely accidentally described. We stress that the 
analyzing power is a very small observable, so that tiny 
improvements of the Hamiltonian can be relevant for a 
proper prediction. This is in line with the rather strong 
dependence on the cutoff, which indicates that order Q 4 
contributions might resolve this puzzle. For the 4N sys- 
tem, a more significant deviation of data a nd predictions 
of the analyzing power has been found by IViviani et al. 
(|200ll ). Interestingly, in this case, chiral 3NFs lead to an 



improved description of the d ata compared t o the stan 



dard phenomenological forces ( Viviani et al. ( 2010l )) 



In summary, the results for few-nucleon systems show 
that Q 3 predictions are in line with the expectations 
based on power counting. Whereas low-energy scatter- 
ing is reasonably described at this order, the results for 
the binding energies indicate that Q A will be required to 
reach satisfactory accuracies. Two-nucleon forces at this 
order are available and h ave an accuracy compa r able t o 
phenomenologica l forc es (jEntem and Machleidtl (|2003l ): 



Epelbaum et all (120051) ). The 3NF s at N^LO h ave been 
completed recently ()Bernard~ ~al\ (j2n08L[201lh V In ad- 



dition, a consistent calculation up to this order also in- 
volves 4N forces, which fortunately do not involve addi- 
tional LECs and are therefore parame ter free. They hav e 
been derived and explored in 4 He (|Epelbauml (|2007l ) ; 



Nogga et all (|2010h ). In this case, for the small cutoffs, 
the contributions seems to be smaller than expected. It 
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remains to be seen whether this is also true for more 
complex systems. 

Next, we return to the discussion of the RG transfor- 
mation started in Section [ITJ because the resolution scale 
dependence also applies to the low-energy couplings in 
3NFs. Therefore, the RG equation in the NN sector, 
Eq. (J7J needs to be augmented by a similar equation for 
3NFs, which we again write schematically as 



— V A (123) = F 123 (V A (ij),V A (U3),A) 



(17) 



For low- momentum interactions V\ ow k, solving the RG 
equation for 3NFs is difficult in practice, because it in- 
volves a compl ete set of scatt e ring s olutions for the three- 
body system (jBogner et all (|2010f) ). This is not feasi- 
ble at this point, but a consistent 3NF evolution can 
be carried out in the SRG approach (as discussed be- 
low). For low- momentum interactions, the chiral EFT 
has been used as a general operator basis of 3NFs, and 
the LECs have been adjusted d irectly to few-nucleo n 
data at lower resolution scales ()Hebeler et Ml (|201lh : 



Nogga et all (|2004)). Such an approach is justified for 



A < 500 MeVj_ because the NN interactions become uni- 



versal ( Bogner et all ((2003)). It is therefore motivated 
that consistent 3NFs should have the same form as the 
ones derived in chiral EFT. Since 3NFs are defined up 
to a finite order, even three-nucleon observables will only 
be approximately independent of the cutoff. It is there- 
fore common to perform calculations for a range of cut- 
off values. The variation of the prediction provides an 
estimate of the theoretical uncertainty due to neglected 
higher-order 3NFs. If observables are calculated for more 
complex systems, the dependence on the cutoff can also 
be due to neglected four- and higher-body interactions. 

The SRG approach provides a powerful scheme to 
evolve 3NFs by differential equations of the general form 
of Eqs. ([7J and (ITT)) . The SRG transformation is an exact 
unitary transformation. Therefore, all NN observables 
are invariant under the transformation. By construction, 
low and high momenta decouple, and observables at low 
mom entum becomes i n sensit ive to high-momentum de- 
tails ( Jurgenson et al. ( 2008h ). As a result, many-body 
calculations converge more rapidly for evolved potentials, 
similar to low-momentum interactions. The SRG evolu- 
tion o f 3NFs has been achieved in a harmonic- oscillator 
basis ( Jurgenson et al. ( 2008 ); Roth et all i20llb ) and re- 
cently in momentum space (jHebelerl ([20121 )). 

The results of the application of SRG-evolved NN 
and 3N interactions to 4 He again shows the quantita- 
tive importance of 3NFs for binding energies of nuclei. 
But it also supports the general belief that four-nucleon 
forces do not contribute significantly, as shown in Fig. 1141 
When the SRG flow is truncated at the two-body level, 
the ground-state energy of 4 He depends significantly on 
the SRG flow parameter A, which plays a similar role 
as the momentum cutoff A. However, the A variation 
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FIG. 14 (Color online) Ground-state energy of 4 He as a func- 
tion of the SRG flow parameter A starting from chiral NN and 
3N interactions at N 3 LO a nd N 2 LO, respectively. For details 
see lJurgenson et al\ (|2009T ) . 



is of the same order as the 3NF contribution. This 
shows how the RG/SRG cutoff variation estimates miss- 
ing parts of the Hamiltonian. When 3NFs are included 
in the SRG evolution, most of the A dependence is re- 
moved. The remaining variation is of the order of 50 kcV 
for A > f.5 fm -1 , indicating that induced 4N forces pro- 
vide a small contribution to the He ground-state en- 
ergy. Note that this estimate is even smaller than ex- 



plicit calcula t ions u sing 4N forces (jDeltuva et all (|2008l ): 
Nogga et all (|201dh ). which result in 200 - 300 keV for 
4 He. 

The small size of 4N forces justifies the exploration 
of larger nuclei and nuclear matter based on chiral NN 
and 3N interactions in the next Section. The detailed 
calculation of 4N forces indicates that the result for 4 He 
might be suppressed because parts of the force cancel for 
these quantum numbers. Eventually, this result needs to 
be confirmed for more complex systems than 4 He. 



V. THREE-NUCLEON FORCES AND MANY-BODY 
SYSTEMS 

Three-body forces are a frontier for understanding 
and predicting strongly interacting many-body systems. 
While the quantitative importance of 3NFs has been well 
established in light nuclei, they arc currently not included 
in most nuclear structure calculations. In this section, we 
discuss the opportunities and challenges this area offers. 
We highlight the importance of 3NFs beyond light nu- 
clei, for neutron-rich systems, and for nucleonic matter 
in astrophysics, with a focus on 3NFs based on chiral 
EFT. Although some of the applications that we discuss 
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FIG. 15 (Color online) Excitation energies in MeV of light 
nuclei, 10 B and 13 C, obtained in the ab-initio No-Core Shell 
Model (NCSM) wi th chiral EFT i ntera ctions (NN to N 3 LO 
and 3N to N 2 LO) (jNavratil et al\ i|2007| y). 



still require an approximative treatment of 3NFs, they 
exhibit new facets and significant contributions of 3NFs. 

As discussed in the previous section, chiral EFT opens 
up a systematic path to investigate many-body forces, 
which has not been possible before. This results from the 
consistency of NN and 3N interactions and the possibility 
to constrain all parameters using only few-nucleon data. 
No new parameters enter for 3N and 4N forces at N 3 LO. 
Moreover, it has been shown that for systems of only 
neutrons, the D and E parts do not contribute because 
of the Pauli principle a nd th e coupling of pions to spin 
(|Hebeler and Schwenkl (f2010h ; iTolos et all (|2008l) ). This 
establishes a forefront connection of the investigation of 
3NFs with the exploration of neutron-rich nuclei at rare 
isotope beam facilities worldwide. 

As expected from the Tjon band in Fig. |H1 3NFs im- 
pact binding energies and therefore also radii. Preci- 
sion techniques for masses and charge radii present new 
challenges for t heory (in the context of 3NFs see, e.g., 



are taken into account. In addition to a repulsive ef- 
fect on the spectra, 3NFs provide important contribu- 
tions to the spin-orbit splitting, reflected in the excita- 
tion energy of the first 3 /2~ state relative to the 1 /2~ 
ground state in 13 C, which probes the splitting of the 
p 3 / 2 — pi/a orbitals. This can also be seen in the 3NF con- 
tributions to the spin-orbit splitting betwee n the pi/ 9 and 
p 3 h ph ase shifts in nucleon- 4 He scattering (jNollett et al. 
(2007)). For 10 B, NN forces alone do not predict the 



correct ground-state spin and parity 3 + , but instead the 
lowest state is found to be 1 + . This is only correc ted by 
some of the phenomenological 3NFs (jPieper (2008)). For 
chiral 3NFs, the correct ordering is predicted. This is also 
needed for the analogous states in medium-mass nuclei 
22 Na and 46 V, which axe N = Z = 8 and N = Z = 20 nu- 
clei with t h ree va lence neutrons and three valence protons 
(jNowackil (|2008l )). Moreover, recent work has demon- 
strated the impact of 3NFs on the structure probed in 
electroweak transitions (see also Section IVII). e.g.. for 
the be t a decay of 14 C used for carbon dating (jHolt et al. 
(|2009h : lMaris et all (|201l h). 



Nuclear lattice simulations were recently used to per- 
form th£_^Jst_abjonMocakulation of the Hoyle state in 
12 C (jEpelbaum et all (|2011l )). which is important for nu- 
cleosynthesis. Due to its alpha-cluster structure, this 
state is challenging for many-body methods. In this ap- 
proach, spacetime is discretized and the nucleons are lo- 
cated on the lattice sites. Their interactions in chiral 
EFT are implemented using auxiliary fields and the low- 
lying states are extracted using a generalized Euclidean 
time projection method. This promising new method al- 
lows to take 3NFs into account without handling large 
interaction matrices. 

The application of RG transformations to evolve nu- 
clear forces to lower resolution leads to greatly en- 
ha nced conve r gence in few- and many-body systems 



(Bog ner et al. 



( 2010h ). Current research focuses on 
extending these methods to 3NFs using the SRG. 
This has been achieve d in a harmonic-oscillator basis 
( Jurgenson et al. ( 2009j ): see Fig. [Hj), with promising re- 
sults in light and medium-mass nuclei in the NC S M an d 
importance - trunca ted NCSM (jjurgenson et all (|201 lh : 
Roth et all (jgOllI) ). and recently in momentum space 



Brodeur et all (|2012l )). In addition, similar to the spin 
dependences observed in few-body scattering, e.g., for 
the analyzing power A y discussed in Fig. (T3J 3NFs play 
an important role for spin-orbit splittings and spin de- 
pendences in nuclei. Both aspects can be clearly seen 
in the spectra of light nuclei, where calculations can be 
performed ab-initio, making these nuclei an interesting 
laboratory to explore nuclear forces. As an example, we 
show two representative spectra in Fig. 1151 calculated in 
the No-Co re Shell Model ( NCSM ) including chiral 3NFs 
at N 2 LO (jNavratil et all (|2007l )). The NCSM is based 
on a large-basis Hamiltonian diagonalization. Without 
3NFs the spectra are generally too compressed (which 
is also found for 23 in Fig. [TH] and for other medium- 
mass nuclei). Clearly, the spectrum improves, when 3NFs stable nuclei (jjanssensl (|2009f ) and references therein) 



(jHebelerl ((2012)). Open questions include understanding 
the cutoff dependence in chiral EFT, whether long-range 
many-body interactions are induced by the SRG, and to 
explore the dependence on the SRG generator. 

Thrce-nucleon forces have also been implemented for 
neutron-rich systems. A frontier in this area is to under- 
stand the sequence of isotopes from proton-rich to the 
limit of neutron-rich nuclei: the neutron dripline. The 
neutron dripline evolves regularly from light to medium- 
mass nuclei e xcept for a s t riking anomaly in the oxy- 
gen isotopes ( Otsuka et all ( 20091) ). where the heaviest 



isotope, 24 0, is doubly-magic and anomalously close to 
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FIG. 16 (Color online) Ground-state energies of the neutron-rich oxygen isotopes relative to 16 0, including experimental 
values of the bound isotopes 16_24 0. Energies obtained from (a) phenomenological forces SDPF-M and USD-B, (b) a G matrix 
interaction and including Fujita-Miyazawa 3NFs due to A excitations, and (c) from low-momentum interactions VJowfc and 
including chiral 3NFs at N 2 LO as well as only due to A excitations. The changes due to 3NFs based on A excitations are 
highlighted by the shaded areas, (d) Schematic illustratio n of a two-valence-neutron interaction generated by 3NFs with a 
nucleon in the 16 O core. For details see lOtsuka et all (|2009h . 
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FIG. 17 Excitation energies of b ound excited state s 
in 23 Q compared with experiment llElekes et all (|2007f ): 
ISchiller et all (|2007l )). The NN-only results are calculated in 
the sd and sdf 7 / 2 P3/2 shells with empirical single-particle ener- 
gies (SPEs). The NN+3N energies are obtained in the same 
spaces, but wi th calculated SPE s including 3NFs at N 2 LO. 
For details see iHolt et all (|201lf) . The dashed lines give the 
one- neutron separation energy S n - 



This oxygen anomaly is not reproduced in shell-model 
calculations derived from microscopic NN forces (see NN 
in Fig. [TC] (b) and (c)), only with phenomenological ad- 
justments (Fig. [T5] (a)). As shown in Fig. [T5] (c), chi- 
ral 3NFs at N 2 LO lead to repulsive contributions to the 
interactions among excess neutrons that change the lo- 
cation of the neutron d ripline from 2 8 O to the exper- 



imentally observed 24 (jOtsuka et all ((2009)). This is 



dominated by the long-range two-pion-exchange part of 
3NFs, as demonstrated in Fig. [1])] (b) and (c) with the 
single- A-excitation Fujita-Miyazawa 3NF (of the type of 



Fig. [TJ. For valence neutrons, the latter contribution is 
repulsive, which can be under stood based on the Pauli 
principle ( Otsuka et al. ( 2009t )). This presents the first 
microscopic explanation of the oxygen anomaly. Since 
the 3NF mechanism is robust and general, these find- 
ings can impact the nucleosynthesis of heavy elements 
in neutron-rich environments. The same 3NF contribu- 
tions have been shown to be key for the calcium iso- 
topes (jGallant et all (|2012f l) an d for valence-proto n in- 
teractions for proton-rich nuclei (H olt et al. (|2012al) ). 

Occupying a position between two neutron-rich, 
doubly-magic isotopes, 22 and 24 0, the spectrum of 
23 in Fig. [IT] provides a unique test for theory, as it 
should reflect the features of both neighbors. In Fig. IT7l 
we observe that 3NF contributions in exten ded valence 
spaces improve the spectrum considerably ( Holt et all 
(|201ll) ). With NN forces, the first excited state is only 
at « 0.5 MeV, well below experiment, simila r to coupled- 
cluste r theory with a N 3 LO NN potential ( Hagen et all 
(120091) ). Important future developments here are to in- 
clude the continuum for loosely bound and unbound 
states ( Michel et al. ( 2010l )). Similar results for the oxy- 
gen and calcium isotopes have recently been obtained 
in coupled-cluster calculations including also the contin- 
uum, where 3NF effec ts are model e d thro ugh adjusted 
two-body interactions (jHagen et al. 

Understanding and predicting the formation and evo- 
lution of shell structure from nuclear forces is another 
key challenge. While the magic numbers N = 2,8, 20 are 
generally well understood, N = 28 is the first standard 
magic number that is not re produced in m i crosco pic the- 
ories with NN forces only (ICaurier et Idl (120051)). In a 
first study for calcium isotopes (jHolt et all (|2012bl) ). it 
was shown that 3NFs are key to explain the N = 28 
magic number, leading to a high 2 + excitation energy 
and a concentrated magnetic dipole transition strength 
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FIG. 18 (Color online) Nuclear matter energy per particle versus Fermi momentum fcp at the Hartree-Fock level (left) and 
including second-order (middle) and third-order particle-particle/hole-hole contribu tions (right), based on evolved N 3 LO NN 
potentials and N 2 LO 3NFs fit to the 3 H binding energy and the 4 He charge radius (jHebeler et al\ l|201lf 0. Theoretical uncer- 
tainties are estimated by the NN (lines) /3N (band) cutoff variations. 



in 48 Ca (|von Neumann-Cosel et aU (|l998l )). 

The calculations of neutron-rich nuclei take into ac- 
count the normal-ordered two-body part of 3NFs, which 
arises from the interactions of two valence neutrons with 
a nucleon in the core (see Fig.Ql)](d)), which is enhanced 
by the number of core nucleons. Moreover, the normal- 
ordered two-body part can be shown to dominate over 
residual three-body interactions ba sed on phase space ar- 
gumen ts for normal Fermi systems (jFriman and Schwenkl 
(l201ll) ). The normal-ordered two-body approximation 
has bee n demonstrated to be effective in coupled-cluster 
theory ([Hagen et all (|2007l) ). with recent ver y promising 



result s for medium-mass closed-shell nuclei (jRoth et al. 

n the context of the shell model, residual 3NFs 



were recently shown to be small, b ut amplified with ne u- 
tron number in neutron-rich nuclei ([Caesar et all ([20121 )). 
In addition, normal-ordering techniques have been used 
to implement the SRG evolution of nuclear Hamiltonians 
directly "in-medium" in the A-body system of interest 
(|Tsukivama et all (|201lh ). 

Recent developments of chiral EFT and RG trans- 
formations for nuclear forces enable controlled calcula- 
tions of matter at nuclear densities. Nuclear matter 
calculations provide an important benchmark for nu- 
clear forces, and are used to constrain calculations of 
heavy nuclei and matter in astrophysics. The RG evo- 
lution to low momenta softens the short-range tensor 
components a n d sho rt-range repulsion of nuclear forces 



(Bog ner et al. 



This leads to contributions in 
the particle-particle channel that are well converged at 
second order in the potential, suggesting that pertur- 
bative approaches can be used in place of the Bethe 



Bruec k ner-Goldstone hole-li ne expansion (jBogner et al. 



also possible to estimate the theoretical uncertainties due 
to neglected many-body forces and from an incomplete 
many-body calculation. The nuclear matter results start- 
ing from chiral EFT interactions are shown in Fig. 1181 
Three-nucleon forces drive saturation, and these are the 
first nuclear forces fit only to A ^ 4 nuclei that predict 
realistic saturation properties. For these developments, 
an improved treatment of 3NFs as density- dep enden t 
two-body interactions has b een key (jHolt et all ([20101 ): 
Hebeler and Schwenkl (|2010l) ). The rapid convergence 
around saturation density in Fig. [18] may justify in part 
the application of in - medium chiral per t urbation the- 
ory ( Lutz et al. (200(3); Kaiser et al. ( 2002 ): Lacour et all 
(|2011[ )). which provides an alternative expansion for nu- 
clear densities. In in-medium chiral perturbation the- 
ory, the inclusion of long-range two-pion-exchange 3NFs 
from A degrees of freedom also improves the descrip- 



tion o f nuclear matter and the convergence (jFritsch et al. 
(2005)). ^n addition, 3NF contributions to the quasipar- 



ticle interactions i n nuclear matter have been explored in 



Holt et all ( 2012cD 



(2005); iHebeler et ail (|201ll )). In this framework, it is 



The nuclear matter results imply that exchange cor- 
relations are tractable, which opens the door to develop 
a universal nuclear energy density functional (UNEDF) 
for global ground-state predictions based on microscopic 
interactions. This is one of the goals of the SciDAC UN- 
EDF/NUCLEI initiatives. Three-nucleon forces play a 
key role in this, including for an improved density ma- 
trix expansion based on chiral EFT interactions (see, 
e.g.. IStoitsov et all (|2010l )) and for systematic studies of 
pairing in nuclei with a non-empirical pairing functional 
t|Lesinski et~al\ (|2012l )). 

For neutron matter, only the long-range two-pion- 
exchange ci and C3 parts of N 2 LO 3NFs contribute 
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TABLE III Symmetry energy E syln obtained from neutron 
matter calculations with N 2 LO 3NFs for different C\ and C3 
couplings and based on RG-evolved N 3 LO NN forces only 
<|Hebeler et al\ l|201Ch ). 



ci [GeV" 1 ] 


c 3 [GeV- 1 ] 


E sym [MeV] 


-0.7 


-2.2 


30.1 


-1.4 


-4.8 


34.4 


NN-onlv CEntem and Machleidt ( 2003)) 


26.5 


NN-onlv CEDelbaum et al. (2005^ 


25.6 



(|Hebeler and Schwenkl ([2010^ ; iTolos et all (|2008l) V This 
has allowed for a detailed study of the theo- 
retical uncertainties o f the neutron matter energy 
(|Hebeler and Schwenkl <|2010ft ). The inclusion of 3NFs 
leads to an energy per particle at saturation density 
E n (p )/N = 16.3 ± 2.2 MeV, where the uncertainty is 
dominated by the uncertainty in the C3 coupling (and 
to a smaller extent by c\ ; see the ci , C3 range in Ta- 
ble HI . Other microscopic calculations lie within this en- 
ergy range. The uncertainty of the prediction is again 
an estimate of the importance of including N 3 LO con- 
tributions for neutron and nuclear matter. Part of the 
N 3 LO 4N forces has been estimated in neutron and nu- 



clear matter (|Fiorilla et all (|2012l) ). and a first complete 
N 3 LO calculation of neutron matter includi ng NN, 3N 



and 4 N forces has recently been carried out (jTews et al. 

The predicted neutron matter energy also pro- 
vides constraints for the symmetry energy (see Table IiTl| 
which demonstrates that the uncertainty in 3NFs domi- 
nates), and predicts the neutron skin thickness of 208 Pb 
to 0.17 ± 0.03 fm, in excellent agreement with a re- 
cent det ermination f r om th e complete electric dipole re- 
sponse (jTamii et These developments are 
complemented by Auxiliary Field Diffusion Monte Carlo 
calculations using a range of phenomenological 3NFs 
( Gandoffi et al. ( 20121) ) and by lattice simulations with 



chiral 3 NFs of dilute neutron matter (jEpelbaum et al. 
(|2009bl n. which can also enable future benchmarks at 



nuclear densities of the perturbative neutron matter cal- 
culations 

These advances have an important impact on astro- 
physics. The microscopic calculations based on chiral 
EFT interactions constrain the properties of neutron- 
rich matter below nuclear densities to a much higher 
degr ee than is refl e cted in current neutron star model- 
ing (Hcb eler et al. Combined with the heaviest 
1.97M neutron star ( Demorest et al. ( 2010f )). the neu- 
tron matter results based on chiral NN and 3N inter- 
actions constrain the radius of a typical 1.4M star to 
R w 10 - 14 km (±15%), as shown in Fig. QU The pre- 
dicted radius range is due, in about equal amounts, to the 
uncertainty in 3N (and higher-body) forces and to the ex- 
trapolation to high densities. The predicted range is also 
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FIG. 19 (Color online) Mass-radius range for neutron stars 
based on chiral EFT NN+3N interacti ons, combined with a 
general extrapolation to high densities (Heb eler et al\ l|2010T ); 
the light blue/grey region includes an update for the 1.97 
neutron star discovered recently). The predicted range is 
consistent with astrophysical modeling of X-ray burst sources 
see, e .g., the red/gray shaded region from the lSteiner et al\ 
20101 ) analysis). For comparison, we also show equations 
of state common ly used in supernova simulations (lines from 
lO'Conneij (|201lD ). 



consistent with astrophysical results obtained from mod- 
eling X-ray burst sources (see, e.g., Steiner et al. ( 2010f) 
in Fig. [T9]) . In addition, the comparisons in Fig. [19] 



demonstrate that the constraints resulting from chiral 
EFT should be included in equations of state used for 
simulations of stellar collapse, black-hole formation, and 
neutron stars. 



VI. THREE-BODY FORCES AND RELATIONS TO 
OTHER PROCESSES 

Because of gauge symmetries, the same expansion is 
used to derive nuclear forces and electroweak operators. 
Therefore, the couplings of three-body forces in an EFT 
determine also electroweak processes. This is an impor- 
tant consistency test and makes such theories very pre- 
dictive. 

A prime example in chiral EFT arc electroweak ax- 
ial currents, where pion couplings contribute both to 
the currents and to nuclear forces. This is already 
seen at leading order: gA determines the axial one- 
body current and the one-pion-exchange potential. Two- 
body currents, also known as meson-exchange currents, 
enter at higher order, just like 3NFs. As shown in 
Fig. [20l the leading axial contributions (at order Q 3 ) 
are du e to long-range on e-pion-exchange and short-range 
parts ( Park et al. (2003)), with the same couplings 03,04 
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FIG. 20 Leading two-body axial currents and the correspond- 
ing 3NF contributions in chiral EFT. Solid (dashed) lines in- 
dicate nucleons (pions), and the wavy line represents the axial 
current. 



and c n of N 2 LO 3NFs (jGardestig and Phillinsl (j2006l) : 
Gazit et aLl (|2009T i). Chiral EFT is essential for this con- 
nection, which can be viewed as the two-body analogue 
of the Goldberger-Treiman relation. 

Two-body current s have also be e n deri v ed for electro- 



magne tic reactions (iPastore et all (12 008): Rollin g et al. 



(2009); 



Rolling et al 



Pastore et al\ (|2009D : IPastore et all (|2011l ): 
(l2011t l). Their application to the few- 
nucleon system has just started, but based on model 
interactions one can expect an interesting sensitivity of 
many electromagne t ic reac t ions to two-b o dy cu rrents and 
3NFs (|Bacca et all (|2009l ): iGolak et all (|2005l )). In this 
colloquium, we focus on recent developments with elec- 
troweak axial currents beyond light nuclei. 

Figure [30] demonstrates the unique constraints chiral 
EFT provides for two-body axial currents and 3NFs. 
This relates the interactions with external probes to the 
strong-interaction dynamics in nuclei. In particular, the 
low-energy coupling C£> that determines the mid-range 
one-pion-exchange 3NF can be determined either from 
the structure of light nuclei (see Section [Vj ; through the 
two-body axial currents that enter weak decays such as 
the 3 H half-life; or from pion production in hadronic 
collisions H This consistency opens up new avenues of 
research for weak interactions and fundamental symme- 
tries. 

Although the importance of two-body currents is 
known from phenomenological studies, for weak pro- 
cesses, chiral currents and the consistency w ith nuclear 



forces h ave only been explo r ed in light nuclei (|Park et al. 



1)20031 ) ; iGazit et all (|2009l) : iRubodera and Rhol ^Ollf TT 
Figure [21] shows the dependence of the 3 H half-life on 
the low-energy coupling cd, which is included both in 
the leading 3NFs and two-body axial currents. With- 
out 3NFs and without two-body currents, the experi- 
mental 3 H half-life is not reproduced. A dependence 
on the different N 3 LO NN potentials is expected, be- 
cause the leading 3NFs and two-body axial currents are 



2 The low-energy coupling cd also enters pion production in NN 
collisions. However, this probes significantly higher momenta, 
because of the produced pion. For nuclear forces, the determi- 
nation in pion production may therefore not be as effective as 
from the low-momentum kinema tics involved in nuclear struc- 
ture (Pandharipandc et al. (200f|)). 
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FIG. 21 (Color online) The ratio (E^) theo /(E^) emp that de- 
termines the 3 H half-life as a function of the low-energy cou- 
pling cd, which relates the leading two-body axial currents 
and 3NFs (see Fig. I20|) . The empirical range is given by the 
horizontal band. Results are shown based on different N 3 LO 
NN potentials and including N 2 LO 3NFs and consistent two- 
body axial currents. For comparison, the result without 3NFs 
and without t wo-body currents ( No MEC, No 3NF) is given. 
For details see lGazit et al\ (|2009t ). 



only order Q 3 . As for 3NFs, the next order two-body 
axial currents are predicted in chiral EFT, without free 
parameters, which enables systematic improvements of 
beta-decay studies and predictions. The chiral EFT 
currents determined from the 3 H half-life have recently 
been applied to the beta decay of th e two-neutron halo 
nucleus 6 He (jVaintraub et all (|2009l )). however using a 
phenomenological potential model not based on chiral 
EFT, where the decay rate is satisfactorily reproduced. 
These theoretical studies are complemented by precision 
measurements (see t he recent result for the 6 He half-life 
(|Rnecht et all ([20121 ))). 

Surprisingly, key aspects of well-known beta decays in 
medium-mass nuclei remain a puzzle. In particular, when 
calculations of Gamow- Teller (GT) transitions of the 
spin-isospin operator g^err 1 * 1 are confronted with exper- 
iment (this is the most significant operator for beta de- 
cays and for electron-capture processes), some degree of 
renormalization, or "quenching" q, of the axial coupling 
9 C A ~ 19a is needed. Compared to the single-nucleon 
value gA — 1.2695(29), the GT term seems to be weaker 
in nuclei. This was first conjectured in studies of beta- 
decay r ates, with a typica l q ~ 0.75 in shell-model calcu- 



lations (jWildenthal et all (|1983I ); iMartinez-Pinedo et al. 



1996)) and other many-body approach e s (|Bender et al. 
20021 ) ; lRodriguez and Martinez-Pinedol (|201fj|) ). In view 



of the significant effect on weak reaction rates, it is no 
surprise that this suppression has been t he target of many 



theore tical works (see the discussion in iVaintraub et al 
(|2009l) ). 
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FIG. 22 (Color online) Nuclear matrix elements M 0vPP for 
neutrino-less double-beta decay of different nuclei. Results 
are shown based on chiral EFT currents at successive orders, 
including one-body currents at orders Q° and Q 2 , and the 
redicted long-range pa rts of two-body currents at order Q 3 
Menendez et all ([201 ll ); for a discussion of the short-range 
contributions, see this reference). For comparison, we also 
show shell-model results (SM09) of IMenendez et all {[200$ ) 
based on phenomenological one-body currents only. 



Recent studies of GT transitions in medium-mass nu- 
clei with chiral EFT currents provide new insi g hts an d 
opportunities to this puzzle ([Menendez et al. 



(2011)). 



Compared to light nuclei, the contributions of chiral two- 
body currents are amplified in medium-mass nuclei be- 
cause of the larger nucleon momenta. Using a normal- 
ordering approximation for two-body currents to create a 
density-dependent operator, it was shown that the lead- 
ing two-body axial currents contribute only to the GT 
operator (up to a small tensor-like correction) and that 
a quenching of low-momentum-transfer GT transitions is 
predicted based on the long-range parts of two-body cur- 
rents. This demonstrates that chiral two-body currents 
naturally contribute to the quenching of GT transitions. 
A reduction of gA in the currents is also expected consid- 
ering chiral 3NFs as density-d ependent two-body interac- 
tions ([Holt et all (|2009l[2010h ). The long-range one-body 
contributions from two-body currents are in part due to 
Delta-hole pairs, but it remains an open problem how 
much of the quenching of gA is due to two-body currents 
and how much due to polarization effects. 

Neutrinoless double-beta decay presents a fundamental 
test of the nature of the neutrino, of le pton number, and 
the ne u trino mass scale and hiera rchy ([Elliott and Vogel 
(|2002T l; lAvignone III et all (|2008f) ). A pivotal input for 
the ongoing and planned experimental searches are the 
nuclear matrix elements that incorporate the structure 
of the parent and daughter nuclei and of the decay 
mechanism. Compared to standard beta decays, neu- 



trinoless double-beta decay prob es different momentum 
transfers Q 100 Me V - m v (jSimkovic et all ([20081 ): 
IMenendez eTaH (|201lh ). Therefore, the rmpact of two- 
body currents is unclear and renormalization effects can 
differ from the suppression of gA in medium-mass nuclei. 
Chiral EFT predicts the momentum-transfer dependence 
of two-body currents, which varies on the order of the 
pion mass due to the one-pion-exchange part in Fig. 1201 
The first calculation of the neutrinoless double-beta de- 
cay operator based on chiral EFT currents at successive 
order is shown in Fig. [2U This demonstrates that the 
contributions from two-body currents are significant and 
should be included in all calculations. It also shows how 
chiral EFT can provide key input and theoretical un- 
certainties for fundamental symmetry tests with nuclei. 
Recently, chiral EFT currents have also been applied to 
calculate the structure factor for spin-dependent weakly 
interacting massive particle (WIMP) sca ttering off nuclei, 
neede d for direct dark matter detection ( Menendez et all 
(l2()12lii. 



VII. OUTLOOK AND FUTURE OPPORTUNITIES 

In this colloquium, we have highlighted the importance 
of three-body forces in nuclear physics and related ar- 
eas. Here, we give an outlook with a focus on future 
opportunities and challenges. Our discussion is guided 
by Fig. [23] which summarizes the leading 3NFs in differ- 
ent EFTs and shows the order in the expansion at which 
they enter. 

In pionless EFT for systems with large two-body scat- 
tering lengths discussed in Section [TTTl three-body forces 
contribute already at leading order because of the Efimov 
effect. If observables are considered at fixed scattering 
lengths, subleading three-body forces are suppressed by 
two orders and enter only at N 2 LO. Some higher-order 
calculations of few-nucleon observables exist but much 
remains to be investigated in this sector. Particularly 
interesting are the application of pionless EFT to halo 
nuclei and low-energy electroweak reactions. Halo nuclei 
are the most promising candidates for observing Efimov 
physics in nuclei, while precise calculations of low-energy 
reactions are relevant for nuclear astrophysics and neu- 
trino physics. In particular, 3NFs play a prominent role 
in two-neutron halo nuclei and larger halo systems. Pio- 
nless EFT also predicts universal three-body correlations 
that can be explored in nuclear reactions in this regime 
and to test the consistency of different theoretical calcu- 
lations (similar to the Tjon line/band). 

In chiral EFT discussed in Sections llV [ IV1 and IV H 
3NFs are suppressed compared to NN interactions. This 
explains the phenomenological success of weaker three- 
body forces of the Fujita-Miyazawa type. As summarized 
in Fig. [23l 3NFs enter at N 2 LO, and their relative contri- 
butions to observables can be understood based on the 



pionless chiral chiral+A 




FIG. 23 Order of 3NF contributions in pionless and chiral EFT and in EFT with explicit A degrees of freedom (chiral+A). 
Open vertices in the last column indicate the differences of the low-energy constants in chiral and chiral+A EFT. 



power counting. Because the operator structure of the 
leading 3NFs is strongly constrained, a global analysis 
of few-body scattering and bound-state data with theo- 
retical uncertainties appears feasible in the framework of 
chiral EFT. This would allow for a determination of the 
long-range Cj couplings in the three-body sector. In addi- 
tion, a consistent determination of two- and three-body 
forces from such an analysis may help to resolve the A y 
puzzle in few-body scattering. 

For applications of chiral EFT interactions to nuclear 
structure, 3NFs play a central role, as discussed for light 
and medium-mass nuclei and for nuclear matter. For 
these many-body calculations, the RG evolution leads to 
greatly improved convergence. A consistent evolution of 
chiral 3NFs has been achieved in a harmonic-oscillator 
basis and recently in momentum space. Important open 
problems are an understanding of the 3NFs induced by 
the SRC and to control higher-body forces, which is nec- 
essary for the desired accuracy in nuclear structure. 

If A(1232) degrees of freedom are included, part of 
the physics contained in the low-energy constants in chi- 
ral EFT is made explicit in lower orders. As a conse- 
quence, a 3NF of the Fujita-Miyazawa type appears al- 
ready at NLO as shown in Fig. [23] Improved convergence 
of the chiral expansion with explicit A degrees of free- 
dom is expected, but a full analysis of few-nucleon data 
remains to be carried out. In addition, a chiral EFT 
with explicit A's would naturally explain why the con- 
tributions from the long-range two-pion-exchange parts 
of 3NFs dominate over the shorter-range parts in appli- 
cations to neutron-rich nuclei and nuclear matter. 

Three- nucleon forces are a frontier in the physics of nu- 
clei that connects the systematic development of nuclear 
forces in chiral EFT with the exploration of neutron-rich 
nuclei at rare isotope beam facilities. The subleading 
3NFs at N 3 LO are predicted in chiral EFT, without free 
parameters, as is the case for N 3 LO 4N forces. In many 
present calculations, the uncertainty of the leading 3NFs 
likely dominates the theoretical uncertainties of the pre- 



dicted observables. The derivation of N 3 LO 3NFs has 
only been completed recently, and no calculation exists 
with N 3 LO 3N or 4N forces beyond few-body systems. 
Therefore, there is a window of opportunity to make key 
discoveries and predictions. In addition to advancing mi- 
croscopic calculations with 3NFs to larger and neutron- 
rich nuclei, an important problem is to understand the 
impact of 3NFs on global nuclear structure predictions, 
e.g., for key regions in the r-process path where system- 
atic theoretical predictions of extreme nuclei, often not 
accessible in the laboratory, are needed. 

Elcctroweak interaction processes are unique probes of 
the physics of nuclei and fundamental symmetries, and 
play a central role in astrophysics. Chiral EFT provides 
a systematic basis for nuclear forces and consistent elec- 
troweak currents, where pion couplings contribute both 
to electroweak currents and to 3NFs. This opens up 
new opportunities for precise nuclear structure calcula- 
tions with theoretical uncertainties that are needed for 
fundamental symmetry tests with beta decays and weak 
transitions, including the key nuclear matrix elements for 
neutrinoless double-beta decay. 

In principle, it is possible to calculate nuclear prop- 
erties directly from the QCD Lagrangian. In Lattice 
QCD, the QCD path integral is evaluated in a dis- 
cretized Euclidean space-time using Monte Carlo sim- 
ulations. This approach is based on a nonperturbative 
formulation of QCD but requires a large numerical ef- 
fort. However, high statistics Lattice QCD simulations 
of two- and three-n ucleon systems are now within reach 
([Beane et all (|2011l )) and the calculation of few-nucleon 
systems appears feasible in the intermediate future. A 
milestone for nuclear forces is the prediction of three- 
neutron properties in a box. This will provide unique 
access to the isospin T = 3/2 component of 3NFs, which 
is not probed in nucleon-deuteron scattering. Moreover, 
Lattice QCD results can also be used to constrain cou- 
plings in chiral 3NFs. A first step in t his direction was 
recently carried out bv lDoi et al. ( 2012 ). 
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